Contour lines of the discrete scale invariant rough surfaces 
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We study the fractal properties of the 2d discrete scale invariant (DSI) rough surfaces. The contour 
lines of these rough surfaces show clear DSI. In the appropriate limit the DSI surfaces converge to 
the scale invariant rough surfaces. The fractal properties of the 2d DSI rough surfaces apart from 
possessing the discrete scale invariance property follow the properties of the contour lines of the 
corresponding scale invariant rough surfaces. We check this hypothesis by calculating numerous 
fractal exponents of the contour lines by using numerical calculations. Apart from calculating the 
known scaling exponents some other new fractal exponents are also calculated. 



I. INTRODUCTION 



Discrete scale invariance (DSI) is a ubiquitous phe- 
nomenon in human made fractals and also in the frac- 
tals of nature. The list of applications of discrete scale 
invariance covers phenomena in a wide range of topics 
such as diffusion in anisotropic quenched random lattices 
[l|, growth processes and rupture 0, non-extensive sta- 
tistical physics econophysics [JTa], soc i°physics @> 
as well as quenched disordered systems Q and cosmic 
lacunarity [8fl, for review of the concept and other appli- 
cations see |9l-[llj . The signature of the presence of DSI 
is complex exponent which manifest itself in log-periodic 
corrections to scaling. In some cases the log-periodic os- 
cillations appear in the time dependence of some physical 
quantities such as the energy release on the approach of 
impending rupture [l2[ and earthquakes [l3| and in some 
other cases the DSI appears in the geometryof the sys- 
tem, the most famous cases being animals [7|, diffusion 
limited aggregation (DLA) 14] and sandpile model on 
Sierpinski gasket [Hj]. Recently we found an other ex- 
ample in surface science, some of the rough surfaces as 
well as their contour lines show clear DSI in their geom- 
etry [l6|. The rough surfaces have been the subject of 
intense studies for many years (l7| . The first predictions 
of the scaling exponents of the contour lines of the scale 
invariant rough surfaces appeared in (l8| . The scaling 
exponents predicted in [18| were confirmed numerically 
in [lj| and the different applications of the findin gs in 
glassy interfaces and turbulence were discussed in [2Cj . 
After the invention of the stochastic Schramm-Loewner 
evolution [2l| which describes scale invariant curves in 
two dimensions a revival appeared in the field and many 
new properties of the contour lines were investigated, for 
example see [2214241. One of the most important find- 
ings of the paper is the dependence of the different 
exponents of the contour lines to the only universal pa- 
rameter in the system, i.e. roughness (Hurst) exponent. 
They showed that for the full scale invariant rough sur- 
faces the fractal dimension of all of the contours d, the 



fractal dimension of one contour Df and the length dis- 
tribution exponent r of the contour lines are just simple 
functions of the roughness exponent. Altough there is no 
theoretical proof for the scaling relations found in [l8| 
they were confirmed in many different numerical studies 

The building block of the functions with DSI is 
Weierstrass-Mandclbrot (WM) function which was stud- 
ied by Berry and Lewis in [25| . It is a random continuous 
non differentiable mono fractal function that one can see 
it in a wide range of physical phenomena such as sedi- 
ment, turbulence (9l.[l0| and propagation and localization 
of waves in fractal media [261 ] . The function is scale in- 
variant just for particular value of the scale parameter 7. 
It was shown in [H, [25[ that in the 7 1 limit the WM 
function converges to the Brownian motion. In the other 
words one can think about WM function as the perturbed 
Brownian motion which the perturbing parameter is the 
7. The same argument is also true for the 2d WM func- 
tion [27], HI] and one can think about the 2d WM function 
as the mono fractal rough surface with the DSI property 
for the specific value of 7. It was shown in [l6j that up 
to the fractal dimension calculations many properties of 
the WM function is the same as the Brownian motion 
counterpart. What we are going to do in this paper is 
confirming this idea for the contour lines of the 2d WM 
function. We will show numerically that many exponents 
of the contour lines of the 2d WM function are the same 
as the exponents of the contour lines of the correspond- 
ing 2d Brownian sheet. In other words we confirm the 
same scaling relations as the paper [l8[ for our discrete 
mono scale invariant contour lines. 

The structure of the paper is as follows: In the next sec- 
tion we will introduce Id and 2d Weierstrass-Mandelbrot 
function and we will discuss the DSI property of the sur- 
faces. The third section is basically for fixing the notation 
and introducing different exponents, most of the materi- 
als of this section were already discussed in 0, . In 
the forth section we will numerically confirm all the pro- 
posed scaling exponents and relations of the third section 
and we will show the universality of those relations. In 
the last section we will summarize our findings and we 
will give some comments about the possible future direc- 
tions. 
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II. DEFINITION OF ID AND 2D 
WEIERSTRASS-MANDELBROT FUNCTION 



The Weierstrass-Mendelbrot (WM) function is a con- 
tinuous non differentiable function. This function has 
discrete scale invariance property which is a weaker kind 
of scale invariance. In other words for the measurable 
quantity O; the relation 0{\x) — p(X)0(x) occurs only 
for a special value of A. The definition of the one dimen- 
sional WM function is 
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(2-D)r 



(1) 



where g(x) is a generic periodic function, 7 > 1 , 1 < 
D < 2 and 4> n is an arbitrary phase [251 ] . This function 
does not change under scaling transformation t — > bt if 
and only if b = 7. In this case WM function has self-afhnc 
scaling law 



W{ji) = 7 2 ~ Z V(<) = i H W(t), 



(2) 



where the parameter H = 2 — D is the Hurst exponent. 
The self affine scaling property of the WM function can 
be also seen in the correlation function of the increment 
AW(t, t) = W(t + t) - W(t) of the process which is 
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where phase <j> n picked up from the uniform random 
distributed numbers. Since the condition C(jt) = 
7 2 ( 2 ~- D )C(r) is true for all periodic functions g(x), there- 
fore the correlation function is also scale invariant. 

The Weierstrass-Mandelbrot function was generalized 
to two dimensions in [23 as follows 



W(x,y) = C 
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where C = (ln 7 /M) 1/2 , B n = (fc o7 ") i -, Vm , n 

1/2 

k cos(9 — a m )j n , r = (x 2 + y 2 ) and 9 = tan -1 (y/x). 
The normalization factor C is chosen so that the cor- 
relation function V (p) be finite in the limits M — >• exo 
and 7—^1. The A m factor can control the isotropy of 
WM surface. We have an isotropic rough surface when 
A m = A. The wave number ko = 77 corresponds to sys- 
tem size L and for each m indices there is one phase a m 
that may be deterministic or stochastic. In this study 
ctm = 77 and (f> m ,n is selected randomly from the inter- 
val [0, 2n] with uniform distribution. 

In numerical methods the two parameters M and n 
should be finite but large enough that W (x, y) does not 
change with increasing these parameters. Fractal dimen- 
sion D and Hurst exponent H = 3 — D are two important 




FIG. 1: (Color online) WM rough surfaces with H = 0.4 
and 7 = 1.5. Light yellow (light Gray) regions correspond to 
larger heights and dark red (dark Gray) regions correspond 
to smaller heights. 



parameters in simulating WM function, for an example 
of 2d WM function see Fig. 1. The surfaces with larger 
H seems smoother than the surfaces with smaller H . 

The correlation function of the two dimensional WM 
function C(p) = I^\W(y + p) — W(v)\ 2 ^j obeys the scal- 
ing law C(-fp) = 7 2< - 3 ~ D - ) C(p) where D is the fractal di- 
mension and 2 < D < 3. This shows that C(p) is also a 
self similar function. 

One of the capabilities of WM function is that some 
of the statistical properties of this function do not have 
any sensitivity to the form of G(x). We use 1 — e lx and 
sin~ 1 (sin(a;)) as G(x) in Eq. @ and we will show that 
the statistical properties of the contour lines of 2d WM 
function depend only to the Hurst exponent of this pro- 



Ill. SCALING LAWS OF THE CONTOUR LINES 



For a given stochastic scale invariant rough surface 
with the height /i(x), a level set /i(x) = ho for differ- 
ent values of ho consists of many closed non-intersecting 
loops. These loops are scale invariant and their size dis- 
tribution is characterized by a few scaling functions and 
scaling exponents. For example the contour line prop- 
erties can be described by the loop correlation function 
Gi(r). The loop correlation function measures the prob- 
ability that two points separated with distance r in the 
plane lie on the same contour. Rotational invariance of 
the contour lines forces Gi(r) to depend only on |r|. This 
function for the contours on the lattice with grid size a 
and in the limit r » a has the scaling behavior 



Gi(r) 
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where xi is the loop correlation exponent. 
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Another measure associated with a given contour loop 
ensemble is G s (r). This probability is called the two 
point correlation function for contour lines with length 
s. Scaling properties of the contours force G s (r) to scale 
with s and r as 



G s (r) 



\r\- n f Gs (r/R), 



(6) 



where to and n are two unknown exponents and 
/gs ( r /R) i s a scaling function. In this function we must 
scale r by the typical diameter R which is called ra- 
dius of gyration. For one loop with N discrete points 
{(xi, yx), (xn, Vn)}, the radius of gyration R is de- 
fined by 



R 2 = 



i N 



Xi - x c f + (yi - y c ) 2 ^j , 



(7) 



where x c = J2i=i x i an d Ve = j? J2i=x Vi are tne cen " 
tral mass coordinates. 

The 2d WM function is a self affine surface. A key 
consequence of this result is that, the contour lines with 
perimeter s and radius R of such surfaces are self-similar. 
When these lines are scale invariant one can determine 
the fractal dimension as the exponent in the perimeter- 
radius relation. The relation between contour length s 
and its radius R is 



(8) 



where Dt is the fractal dimension, R is defined by Eq. 
([7]) and s is measured with a ruler of length a. It is a right 
time to mention that the Df is the fractal dimension of 
one contour and it is different from the fractal dimension 
of all the level set d. The latter one can be derived from 
the fractal dimention of the rough surface which is 3 — H 
and the intersection rule of the Mandelbrot. Since the 
level set, the set of the contour lines, is the intersection 
of the rough surface with a two dimensional surface one 
can derive the following relation 



d=2-H. 



(9) 



For a given self similar loop ensemble one can define 
the probability distribution of contour lengths P(s). This 
function is a measure for the loops with length s and 
follows the power law 



P(s) 



(10) 



where r is a scaling exponent. The above power law func- 
tions (Eqs. ©, © and (fl"U)) ), introduce the geometrical 
exponents xi, Df and r. These exponents depend only 
on the roughness exponent H (l8| . 

One can derive mean loop length < s > and the prob- 
ability distribution of length P(s) from another quantity 
n(i?, s). This function is the probability distribution of 
contours with length (s, s + ds) and radius (R, R + dR). 



Scale invariance of the contours forces h(R, s) to scale 
with s and R as 



(11) 



In the above relation y is related to the fractal dimension 
Df and the length distribution exponent r. It is obvious 
that 



<s>(R) 



P(s) 
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J °° dsh(R, s) ' 



dRh(R, s). 



(12) 
(13) 



The above equations and Eq. (jlll) yield P(s) ~ 
s~ v+1 / D f j which we can equate this with Eq. (|10p. and 
obtain y = r + 1/Df. Notice that, one can derive scaling 
law © from Eqs. CO]) and dHJ. 

Following [18] one can determine the relation between 
all geometrical exponents. In order to find the first scal- 
ing relation, one can use from the scale invariance of the 
rough surfaces, i.e. h(bx) = b H h(x) under x — > bx, 
and invariance of n(R) under rescaling, i.e. n(R/b) = 
b 3 ~ H n(R), then 



h(R) - R 



-3+H 



(14) 



On the other hand, the number of loops with radius in 
(R, R + dR) can be derived from h(R) = J °° dsn(R, s), 
where it obeys the scaling law 



n(R) ~ R-V+W-W. 



(15) 



Equations (fT4|) and (|15|) lead to the first scaling relation 



D f (r-1)=2-H. 



(16) 



This relation is called hyperscaling relation; it is derived 
for the scale invariant rough surfaces. The above argu- 
ment is also true for the discrete scale invariant surfaces 
such as 2d WM functions. For the DSI rough surfaces 
such as Eq. (j4|), condition h(bx) — b H h(x) is true just 
for a special value of b = 7. 

For a self affine random field, there is another scal- 
ing relation that connect the three geometrical exponents 
Df , t and x\. This relation is derived from a sum rule 
[18| . Integrating (j6]) over all lattice points gives the total 
number of points on the loop 



s = J d 2 rG s (r) 



(2-n)/Df+m 



(17) 



The above equality gives one relation between the expo- 
nents n and m as 



2-n = D f {l - m). 



(18) 



The integration of G s (r) over loop probability distribu- 
tion function P(s) = sP(s) gives the total loop correla- 
tion function Gi(r) 

G t (r) = [ dsP(s)G s (r) ~ r~ w (r^ ) m+2 ~ T , (19) 
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where after using Eqs. (J5)) and (Ti~9t gives another relation 
between m and n as 



D f (m + 2 - r) 



-2xi. 



(20) 



Elimination of m and n from equations (|18l) and (|20[) 
leads to the second scaling relation 



D/(3 - t) = 2 - 2.t ; . 



(21) 



If we combine the two scaling relations (jTHJ) and (j2"Tj) we 
can find two exponents Df and r as a function of the 
Hurst exponent H and the loop correlation exponent xi , 



and 



D i 
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H 
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In the next section we will numerically calculate all 
the introduced exponents related to the contour lines of 
DSI rough surfaces and we will show that the two scal- 
ing relations (|T6|) and (|2~T1) are valid for generic 2d-WM 
function. 



IV. NUMERICAL RESULTS 



To simulate the 2d WM function on the discrete lattice 
we used lattices on square grid with sizes 1024 and 2048. 
In equation (j4j we set M = 20 and —n m i n — n max = 
50. After setting 7= 1.5 we simulated the WM rough 
surfaces with the Hurst parameter in the range of H = 
{0.0,0.1,0.2,0.3,0.4,0.5,0.6} for two periodic functions 
G(x) = (1 — e lx ) and G(x) — sin (sin(;r)) . In each 
case, all numerical tests have been done using from 250 
realizations. The contouring algorithm followed from [23] 
to generate contour lines of rough surfaces. In Fig. 2 one 
can see that the number of loops increases by decreasing 
roughness exponent. 

We would like to measure the exponents d, Df, r and 
xi with a large enough loop ensemble for the mentioned 
values of the Hurst exponents. Before starting numerical 
computation of the geometrical exponents, it is impor- 
tant to see the DSI property in the contour lines of such 
discrete scale invariant rough surfaces. At the beginning 
we will show that the contour lines of DSI rough sur- 
faces are discrete scale invariant. Existence of DSI in the 
loop ensemble can be studied by the Lomb normalized 
priodogram test. 



by using periodic function G(x) with infinite number of 
discrete frequencies u) n = 7™. Each frequency contributes 
in (j4j with weight j~ Hn . To detect these frequencies one 
can use the Lomb normalized periodogram. 

For a given time series Xi = X(ti) with N data points, 
this kind of spectral analysis fitting these data points 
to a sine function of varying frequencies by least square 
method (29j . Lomb analysis measures spectral power 
P(uj) as a function of angular frequency uj as 



1 f \E j (X j -X)coau}(t j -8)] 2 
2a 2 \ J2j^s 2 cj(tj -6) 



V sin u(tj - S) 



(24) 



where X= i£f X h a 2 



X) 2 and 6 is 



1 r* i i i i /c\ c\ sin2cut 7 - 

denned by the relation tan(2w<3) = ^ cos 2u t- • 

We used from Lomb test to show that the contour 
lines of the DSI rough surfaces are discrete scale in- 
variant. For a closed loop (see Fig. 3), with TV data 

points {(xi,j/i), (xn, Hn)}-, the distance from the cen- 

1/2 



ter of mass; i.e. 



r(i) = ((a 



,r + {yi - vcY 



where x c = jj X)i=i x i an d Vc — jf J2%=i Vi are tne cen " 
ter of mass coordinates; is a good quantity to detect 
the discrete frequencies. Using Lomb test we measured 
lexp = (^i+iA^i) = l-55±0. 04 where Wj+i and are two 
frequencies that P(oS) has sharp picks, results are shown 
in Fig. 4. This result is in a good agreement with the 
theoretical value of 7 = 1.5. 



B. Loop correlation exponents 

In this subsection we will numerically find the most 
fundamental exponents of the contour lines of the rough 
surfaces, i.e. xi and n, by analyzing the functions G s (r) 
and Gi(r). 



1. Scaling exponents in Gi(r) 

To find the correlation function from a given loop 
ensemble we follow the algorithm described in [18|]. Us- 
ing scaling properties of loop correlation function we can 
measure correlation exponent x\. For the large class 
of rough surfaces xi = |, this scaling exponent is in- 
dependent of H [l8{. Our numerical tests show that 
xi = 0.5 ± 0.05 for all the WM rough surfaces and it 
is not a function of H or G(x) as shown in Fig. 5. 



A. The Lomb normalized periodogram 



2. Scaling exponents in G s (r) 



The 2d Weierstrass Mandelbrot function (Q} is a dis- 
crete scale invariant function which can be constructed 



The scaling behavior of the two point correlation func- 
tion of the contours with length s is defined in Eq. ([5]). 



FIG. 2: (Color online) Small parts of contour lines of 2d WM rough surface with G(x) = (1 — e lx ) and lattice size 1024 2 . Top 
left (H = 0.0), right (H = 0.2), and bottom left (H = 0.4), right {H = 0.6). 




FIG. 4: (Color online) Lomb test on the contour line of 2D 
FIG. 3: (Color online) Largest contour line of 2d WM rough WM surfaces wit h D = 2.5 and 7 = 1.5. Lomb test on Id 
surfaces with N = 6 x 10 3 , D = 2.5 and 7 = 1.5. Distance WM time series with y = 4 0) is shown in the inset sharp 
r(i) for each point i from the center of mass is shown in the picks show that the contour lines of WM sur f aces are DSL 
inset. 



The exponent n can be measured from the power law de- 
pendence of G s (r) with respect to r for a fixed value of 
the perimeter s, see Fig. 6. For this purpose the contour 
loops with the perimeter in the range (s, as) (a > 1) 



are considered and G s (r) is measured as a function of r 
for the different values of s, see Table 1. There is nei- 
ther theoretical prediction for this exponent and nor any 
numerical estimates. Our numerical calculation shows 
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FIG. 5: (Color online) The loop correlation function for the 
DSI rough surfaces with H — 0.4, 0.5 and 0.6 (from top to 
bottom) for G(x) = 1 — e !a: . The loop correlation function for 
G(x) = sin _1 (sin(a;)) is shown in the inset. The plots were 
derived for the system size L — 2048 and 7 = 1.5. 



that the exponent n closely follows the relation i4j— ■ , see 
Fig. 7 which means that the exponent m is very small 
compared to the exponent n. Our numerical calculation 
shows that the same values can also be derived for the 
scaling mono fractal rough surfaces. 




FIG. 6: (Color online) The scaling behaviour of G a (r) as a 
function of r in the DSI rough surfaces for G(x) = 1 — e lx 
with L = 2048, H = 0.5 and 7 = 1.5. 



C. Fractal dimensions 

To study the fractal properties of contour lines, the 
fractal dimension of loops D / and the fractal dimension 
of all the contours d, in DSI rough surfaces one can use 
from three related approaches. In the first method one 
can use from self-similar properties of contours to mea- 
sure fractal dimensions. The second method is using cu- 
mulative distribution of area and the last one is the Zipf 's 
law. We used these three methods to find the two geo- 
metrical exponents Df and d. 



FIG. 7: (Color online) The exponent n with respect to the 
roughness exponent H for G(x) = l — e tx . The slop of the best 
fit line is 0.45 ± 0.03. The smooth line is the line n 



l+H 
2 



H 


G(x) = l-e rx 


G(x) = sin (sin(a;)) 


0.0 


0.46 ± 0.04 


0.47 ±0.03 


0.1 


0.53 ±0.03 


0.52 ±0.03 


0.2 


0.56 ±0.04 


0.56 ±0.02 


0.3 


0.59 ±0.03 


0.59 ±0.04 


0.4 


0.64 ±0.03 


0.64 ±0.03 


0.5 


0.69 ±0.03 


0.69 ±0.04 


0.6 


0.74 ±0.03 


0.75 ±0.03 



TABLE I: The numerical values of the scaling exponent n 
derived from the two point correlation function G s (r). 



In order to find the fractal dimension of the contours 
Df we used directly from the scaling relation between 
loop perimeter s and the radius of gyration R accord- 
ing to Eq. (|SJ). Scaling behavior of the averaged loop 
length < s > with respect to R in log-log plot is shown 
in Fig. 8. The geometrical exponent Df is measured 
using the least-square linear fit in the scaling regime 
(10 < R < 100). 

The fractal dimension of_all contours d can be found 
by box counting method [3l|. This method is based on 
the self similar properties of all contours in the level set. 
For a self similar object there is a scaling law between 
the number of boxes N(L) with size L needed to cover 
the fractal object 



N(L) - L~ 



(25) 



a. Length-radius scaling relation: 



As shown in Fig. 9, we find excellent agreement between a 
measured values of Df and d with theoretical predictions. 

b. Cumulative distribution of areas: 

The cumulative distribution of the loop area is a good 
measure to describe the fractal dimension of the set of 
all contour lines d 23]. For a self affine random field 
the number of contours with area greater than A has the 
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FIG. 8: (Color online) The scaling relation between the av- 
eraged loop length < s > and the radius of gyration R for 
contour lines of the DSI rough surfaces with G(x) = 1 — e lx . 
The same scaling behavior for G(x) = sin^^sinfx)) is shown 
in the inset. 
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FIG. 9: (Color online) The fractal dimension of all contours 
d and the geometrical exponent Df as a function of H for the 
DSI surfaces with G(x) = 1 — e rx . These scaling exponents 
for G(x) = sin -1 (sm(a;)) are shown in the inset. The solid 
lines correspond to the theoretical predictions. 

scaling form 

N>(A) ~ A- d/2 , (26) 

where d = 2 — H . An example of the scaling behaviour in 
cumulative distribution of area is shown in Fig. 10. Here, 
we plot N > (A) as a function of A for H = 0.4 in log-log 
plot, to show how they follow Eq. (f26|) . In Table ITT1 we 
report the scaling exponent d/2 for different values of H. 

c. Zipf's laws 

Zipf's law is one of the most prominent examples of 
power laws in nature 0, I3(il - l33l | . If we assign ranks to 
a measurable quantity Ai n according to their size, there 
is an interesting power law relation in the size-rank dis- 
tribution which is called by Mandelbrot the Zipf's law. 
In this method the greatest one measured quantity has 
rank 1 (Mi), smaller rank 2 (M2) and so on. Zipf's law 
states that the ranked size M. n falls off as 

M n ~ n-i, (27) 
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A 

FIG. 10: (Color online) Cumulative number of loops whose 
area is greater than A for the DSI rough surfaces with G(x) = 
1 — e' x and G(x) — sin~ 1 (sin(a:)) with H — 0.4, system size 
L = 2048 and 7 = 1.5. 
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Theory 


G(x) =1-6" 


G(x) = sin 1 (sin(a;)) 


0.0 


1.00 


1.00 ±0.01 


1.00 ±0.01 


0.1 


0.95 


0.96 ± 0.01 


0.96 ±0.01 


0.2 


0.90 


0.91 ±0.01 


0.90 ±0.01 


0.3 


0.85 


0.86 ± 0.01 


0.86 ±0.01 


0.4 


0.80 


0.81 ±0.01 


0.81 ±0.01 


0.5 


0.75 


0.76 ± 0.01 


0.77 ±0.02 


0.6 


0.70 


0.71 ±0.01 


0.71 ±0.01 



TABLE II: The best fit values of the exponent d/2 derived 
from the scaling laws of cumulative distribution N > (A). 

where 7 is a scaling exponent. The scaling laws (|27|) for 
■M-n = {'n, A n , R n } with l n , A n and R n are in order, 
the ranked perimeter, ranked area and ranked radius of 
gyration, arc 

l n ~n- D f/ d , A n ~ n - 2 ' d , R n ~n- l l d , (28) 

where Df is the fractal dimension of contours and d is 
the fractal dimension of all loops. We calculated the ex- 
ponents in Eq. (|28|) for 2D WM rough surfaces using 
numerical tests. Averaged values over different realiza- 
tions for each exponents in Eq. (|28]) as a function of H 
is presented in Figs. 11. 

D. Length distribution exponent 

The loop exponent t can be extracted from P(s) which 
follows Eq. (fTUjl . This probability distribution function 
at large s has the statistical noise. In order to minimize 
these noises, we consider the cumulative distribution of 
loop length P>(s) = f°° P(s')ds' with the scaling form 
P>( s ) ~ s~ T+1 - The length distribution exponent r can 
be measured numerically from the power-law scaling in 
loop length (10 < s < 1000). The scaling behavior of 
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FIG. 11: (Color online) The numerical estimation for the scal- 
ing exponents: top 1/d, middle Df /d and bottom 2/d; for the 
WM rough surfaces with G(x) = 1 — e lx and the system size 
L = 2048. The numerical computation of each exponents for 
the rough surfaces with G(x) — sin _1 (sin(x)) are shown in the 
insets. The solid lines show the theoretical relation between 
1/d, Df/d and 2/d with the Hurst exponent H. 



the cumulative distribution -P>(s) for H = 0.5 is plot- 
ted in Fig. 12 for different periodic functions G(x). The 
deviation from scaling line at large s comes from finite 
size effects. Numerical results for the exponent 1 — r are 
summarized in Table ITTTl 



V. CONCLUSION 

In this paper we first introduced a wide range of rough 
surfaces with the discrete scale invariance property ((H). 
The surfaces in the limit 7 — > 1 converge to the Gaussian 
mono fractal rough surfaces. The contour lines of the 
introduced curves show clear DSI property. The frac- 
tal dimension of the contours follow the behavior of the 
fixed point, i.e. Brownian sheet. The numerical calcu- 



10 1 10 2 10 3 10 4 10 5 

s 

FIG. 12: (Color online) Cumulative loop-size distribution, 
for DSI rough surfaces with H — 0.5 and system size L = 
2048. Note that this plot shows that the scaling exponent r 
is independent of G(x). 



H 


Theory 


G{x) = l-e ix 


G(x) — sin 1 (sin(x-)) 


0.0 


-1.33 


-1.33 ±0.01 


-1.34 ±0.01 


0.1 


-1.31 


-1.31 ± 0.01 


-1.32 ±0.01 


0.2 


-1.29 


-1.29 ± 0.01 


-1.29 ±0.01 


0.3 


-1.26 


-1.27 ±0.01 


-1.26 ±0.01 


0.4 


-1.23 


-1.24 ±0.01 


-1.24 ±0.01 


0.5 


-1.20 


-1.21 ± 0.01 


-1.21 ±0.01 


0.6 


-1.17 


-1.18 ± 0.01 


-1.18 ±0.01 



TABLE III: The best fit values of the exponent 1 — r derived 
from the scaling laws of cumulative distribution P>(s). 

lation shows that the relation xi = \ is superuniversal 
which means that it is independent of G{x), H and 7. 
Insteed the exponent n changes linearly with respect to 
the roughness exponent H . We also checked numerically 
the consistency scaling exponent relations, i.e. (|16[) and 
(f2Tj) . There are different methods to define WM rough 
surfaces in two dimensions, for a different example see 
[2^ |. we strongly believe that the same scaling relations 
also hold for these surfaces. Although the relations that 
we introduced in section 3 are valid for almost all the 
mono fractal surfaces the careful studies of the multi frac- 
tal surfaces shows huge discrepancies [3~H - 

Our calculation shows that it is possible to study the 
DSI version of the contour lines of the scale invariant 
rough surfaces. The DSI property does not change the 
fractal dimension of the curves. It seems that this is 
a very general property [lq and possibly can be seen 
in a wide range of statistical models. In other words it 
should be possible to change the usual statistical models 
in a particular way to have DSI property with the similar 
fractal properties as the ordinary statistical mechanics 
models. 

Up to know the only difference of our models with the 
scale invariant rough surfaces is the DSI property. It is 
quite natural to try to find some other quantities that can 
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distinguish the continuous scale invariant rough surfaces 
from the discrete scale invariant rough surfaces. Since 
most of the well-known statistical models at the critical 
point show universality and in some cases conformal in- 
variance it seems that we need a better understanding of 
discrete scale invariant field theories. 
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